Weak coupling BCS suerconductivity in the Electron-Doped Cuprate Superconductors 
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We use in-plane tunneling spectroscopy to study the temperature dependence of the local su- 
perconducting gap A(T) in electron-doped copper oxides with various T c 's and Ce-doping con- 
centrations. We show that the temperature dependence of A(T) follows the expectation of the 
Bardeen-Cooper-Schrieffer (BCS) theory of superconductivity, where A(0)/fesTc ~ 1.72 ± 0.15 and 
A(0) is the average superconducting gap across the Fermi surface, for all the doping levels inves- 
tigated. These results suggest that the electron-doped superconducting copper oxides are weak 
coupling BCS superconductors. 

PACS numbers: 74.50.+r, 74.72.Jt, 74.45. +c 



I. INTRODUCTION 



The physics of conventional superconductors can be 
well understood by the Bardeen-Cooper-Schrieffer (BCS) 
theory of superconductivity. Within the BCS model, 
the superconducting gap A(T) is weakly temperature de- 
pendent at low temperatures but closes rapidly to zero 
near T c . In the weak coupling limit, A(0) /ksT c is 1.76 
for an isotropic gap and becomes slightly smaller for an 
anisotropic gap, where A(0) is the zero-temperature gap 
averaged over the entire Fermi surface. Since there is no 
generally accepted microscopic theory for high-transition 
temperature (high-T c ) copper oxides, it would be inter- 
esting to see if the BCS theory can under certain condi- 
tions describe the physics of some high-T c cuprates. For 
hole-doped (p-type) materials, a pseudogap appears at 
the antinodal region and may compete with the supercon- 
ducting gap on the Fermi surface [3,13. Electron-doped 
(n-type) copper oxides, on the other hand, have a much 
weaker pseudogap effect and thus provide a good oppor- 
tunity to investigate the superconducting gap without 
the influence of "Fermi arc", "Nodal metal" or "Pseu- 
dogap" [E S H, S Q • Using angle resolved photoemis- 
sion spectroscopy and transport measurements, previous 
work have found that the Fermi surfaces in the n-type 
cuprates have two-band characteristics @, 0, H, @, E3 and 
the superconductin g gap has the non-monotonic d-wave 



pairing symmetry 
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}. Unfortunately, the super- 
conducting gap of n-type cuprates is relatively small and 
determination of its exact value over a large doping range 
is an experimental challenge [3, EU EH [HI, 03 Ell ■ Pre- 
vious estimates suggest that the superconducting pair- 
ing strength of the n-type cuprates is close to a weak- 
coupling regime in the optimally doped and overdoped 
region [18l.ll9j|. Furthermore, tunneling data suggest that 



TABLE I: Main superconducting phases studied in this work. 



Label 



Annealing T e (K) 



Formula 



plcco-un21 underdoped 21 

plcco-un24 underdoped 24 

ncco-op25 optimally doped 25 

plcco-ovl7 overdoped 17 

plcco-ovl3 overdoped 13 



Pro.ssLaCeo.i2Cu04- y 
Pro.ssLaCeoA2CuO4.-y 

Nd1.s5Ceo.mCuO4.-y 
Pro.s5LaCeo.i5Cu04- y 
Pro.si,LaCeo.i5,Cu04-y 
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the superconducting gap increases monotonically with 
decreasing doping levels even in the underdoped regime 
[rr| . Although much is known about these electron- 
doped cuprates, there have been no systematic study of 
the superconducting gap as a functon of electron-doping. 

In this paper, we present the point contact spectra mea- 
sured on n-type cuprates Nd2- x Ce x CuOi-y (NCCO) 
and Pri- x LaCe x CuOi- y (PLCCO) over a wide doping 
range. We find that the temperature dependence of the 
average gap follows the BCS predictions with a universal 
weak coupling ratio for all doping levels investigated. 



II. EXPERIMENTAL 

High-quality single crystals of NCCO and PLCCO 
were grown by the traveling-solvent floating-zone tech- 
nique. As grown, the crystals are not superconducting. 
By annealing the samples at different temperatures in 
pure Ar or vacuum, bulk superconductivity with differ- 
ent T c 's can be obtained [20J. The detailed information of 
the superconducting phases studied in this work are pre- 
sented in Table |TJ The in-plane point contact junctions 
were made by approaching the Pt/ Ir alloy or Au tips to- 
wards the (100) and/or (110) surfaces of the single-crystal 
samples (As discussed in this paper, there is no obvious 
difference between these two directions). The tip's prepa- 
ration and the details of the experimental setup were de- 
scribed elsewhere [2l[. In order to obtain high-quality 
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FIG. 1: Point-contact spectra measured at 2K and well above 
T c for various doping levels. The conductance drop has been 
subtracted from the spectra above T c as explained in the text. 
Those two vertical grey lines are added as references to iden- 
tify the positions of the coherence peaks. 



junctions with good reproducibility, the samples were 
carefully processed by nonaqueous chemical etching be- 
fore being mounted on the point contact device [l(| ■ For 
each superconducting phase, we repeated measurements 
many times at different locations on the sample surface 
and obtained the local superconducting gap A and tran- 
sition temperature T c which vary slightly around the bulk 
values. 



III. RESULTS AND DISCUSSIONS 

Fig. [1] shows the point contact spectra measured at 
various doping levels. These spectra are highly repro- 
ducible on the sample surface. In each run of the mea- 
surement, the spectra were recorded at various tempera- 
tures between 2 K and T c with increments of 1 K (refer 
to Fig.[3Ia)-(c)). The data measured at T = 2 K and well 



above T c are presented in Fig. [T] for clarity. On the low 
temperature spectra, two coherence peaks are accompa- 
nied by low-energy depression of the quasiparticle density 
of states. The conductance within the gap voltage does 
not go to zero because the junctions are not tunnel junc- 
tions but ballistic point contact junctions with a finite 
tunneling barrier [22l . |23| . One puzzling aspect of the 
data is that no any zero bias conductance peak (ZBCP) 
is found along the nodal direction for all doping levels, in- 
consistent with the expectation of non-monotonic d-wave 
pairing symmetry for n-type cuprates [ill EH El> HH • In- 
stead, the spectra show almost identical shape with two 
distinct coherence peaks. This is difficult to understand 
within the current d-wave theory. This may be caused by 
the micro roughness (or corrugation) of the nodal surface 
arising from the inherent crystal-lattice structure and the 
axis dependent strength of chemical bonding. In this 
case, the dominated incident current is actually injected 
along antinodal direction of the sample, and hence the 
identical spectral shape for both nodal and antinodal di- 
rections can be easily understood. 

It is well known that when temperature rises across the 
superconducting transition, the junction becomes normal 
and the background conductance drops notably due to a 
finite normal-state resistance of the sample (R s ) in series 
with the junction (24[. In this work, R s varies from 
to 2 ohm depending on the samples' property and the 
configuration of the electrodes and the point contacts. 
It can be subtracted simply by replacing dV/dl and V 
by dV/dl — R s and V — IR S , respectively. All spectra 
(above T c ) shown in Fig. [T] have been treated in this way. 
Accordingly, we can readily construct a universal back- 
ground and normalize the spectra below T c which were 
then compared with theoretical models. In addition to 
the rough estimation of 2A from the peak-to-peak dis- 
tance of the spectra (it is well known that such method 
always overestimates the gap value when the measured 
spectrum is not ideal, in this work such overestimate can 
exceed 1 meV) , we used the BTK theory [22j to derive A 
more accurately [H, H|| . In this model, two parameters 
are introduced to describe the effective potential barrier 
(Z) and the superconducting energy gap (A). As a sup- 
plement, the quasiparticle energy E is replaced by E+iT, 
where T is the broadening parameter characterizing the 
finite lifetime of the quasiparticles [U HI] . For the ex- 
tended anisotropic BTK model [II], another parameter 
of a was introduced to distinguish different tunneling di- 
rections. In this work, a is set to corresponding to the 
antinodal direction as discussed above. 

In Fig. Ufa) and (b), we present the comparison be- 
tween the experimental data and the theoretical cal- 
culations for both isotopic and anisotropic BTK mod- 
els. Fig. [2fc) illustrates the gap function of the non- 
monotonic d-wave pairing symmetry expressed explicitly 
by A ams = A o [1.43cos(20) - 0.43cos(6^)] [H, EH- In 
some cases, the nodal region (the grey segments) has neg- 
ligible contribution compared to the antinodal region (the 
black segments) because of the weaker spectral weight [1] 
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FIG. 2: Comparison of experimental data and theoretical cal- 
culations with the isotropic BTK model (black lines) and the 
anisotropic one (grey lines) for two samples: (a) plcco-un21, 
(b) ncco-op25. (c) The schematic diagram of the anisotropic 
gap (or non-monotonic d-wave gap), in which the black frag- 
ments indicate the antinodal region and the grey fragments 
indicate the nodal region, (d) A(averaged gap)~ T relations 
determined with two different models, solid lines are guides 
to eyes. 



or fewer charge carriers d, [nj possessed by the nodal 
Fermi pockets. As the first order approximation, we only 
consider the contribution from the antinodal region and 
describe the weakly angle dependent gap as a constant 
value, which is equal to the isotropic BTK model. As 
shown in Fig. HJa), both the isotropic BTK model and 
the non-monotonic d-wave one fit the data very well, pos- 
sibly due to the broadening effect (a finite T value). If 
we look at the spectra shown in Fig. &b) which has a 
smaller broadening effect (T/A w 0.28) than that in 
Fig.^a) (r/A rs 0.44), the non-monotonic d-wave model 
fits the data better at higher energy outside the coher- 
ence peaks while the isotropic model is more favorable 
at lower energy. This is consistent with the above dis- 
cussions, i.e., the contribution from the antinodal Fermi 
pockets is dominant in this spectrum and the constant 
gap is a good approximation. In Fig. [2Id), we present 
the temperature dependence of the averaged supercon- 
ducting gap determined by both models. We find that 
the averaged gap determined by the isotropic model has 
a very small uncertainty below 10% for all the doping 
levels. 

Fig. [3] illustrates the details of our data analysis. 
Fig. C2a)-(c) show the temperature dependence of the 
normalized spectra of an underdoped sample, the op- 
titimally doped one and an overdoped one, respectively. 
These data are consistent with the calculations based on 
the isotropic BTK model (denoted by the solid lines). 
As shown in Fig. OHd)-(i), the value of the fitted barrier 
strength (Z) lies between 1.0 ~ 1.5 and Y value changes 
in a range from 0.7 meV to 1.5 meV without an obvious 





Underdoped: T c =19.5, z)(0)=3.10 
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FIG. 3: (Color online) Temperature dependence of the nor- 
malized spectra (balck dots) and theoretical calculations with 
the isotropic BTK model (solid lines) for: (a) plcco-un21, (b) 
ncco-op25, and (c) plcco-ovl7, all curves except the lowest one 
are shifted upwards for clarity; (d)-(f) Fitting parameter of Z 
corresponding to the data shown in (a)-(c) respectively; (g)- 
(i) Fitting parameter of T corresponding to the data shown 
in (a)-(c) respectively; (j) The temperature dependence of the 
superconducting gap in a reduced scale for various doping lev- 
els. 



dependence on the doping levels. The independence of Z 
on temperature indicates high stability of the junctions. 
For the optimally doped and overdoped samples, T is also 
almost independent on temperature and increases with 
junction resistance similar to that of conventional BCS 
superconductors such as Nb [3(| and Zn [3l|. However, 
the T-value of the underdoped sample decreases continu- 
ously when temperature increases very close to T c . The 
temperature dependence of T may be closely related to 
the inhomogeneity, impurities, disorders, and scattering 
mechanism in this system, which need to be clarified by 
future experiments. In fact, it has been demonstrated 
that the T value is not an obstacle to derive the gap value 
for a wide temperature scope [To| . For example, when T 
just exceeds A (r/A « 1.1), the uncertainty of A is still 
below 30%. However, there is another factor making the 
fitting procedure more difficult, which was called as "crit- 
ical current effect" [j^, because it will distort the shape 
of the spectra around T c . Therefore, the upper temper- 
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FIG. 4: Relationship between the zero-temperature averaged 
gap and transition temperature. The grey symbols indicate 
previous reports from ARPES 13], Raman [12], point contact 
PUS an d grain-boundary [lj], respectively. The open cir- 
cles denote the results of fitting to the non-monotonic d-wave 
model. The crosses indicate the data of Nb measured by the 
same experimental apparatus used in this work. The data of 
some conventional superconductors are also presented as open 
stars for comparison [29| . 

ature limit of our analysis is mainly determined by such 
effect. Accordingly, we gave up to fit the spectra for the 
temperatures very close to T c , as presented in Fig.[3fa)- 
(c). ^ 

It is generally accepted that a practical point contact 
junction often includes many channels or real point con- 
tacts, so it is difficult to estimate the contact area simply 
according to the Sharvin formula. Alternatively, we have 
taken a prudential way to keep our measurements away 
from artificial errors. As elaborated in previous reference 
[HI], after the metal tip reaches the sample surface, the 
barrier layer is abraded at first and its thickness decreases 
slightly with the increasing pressure. Consequently, 
the measured spectrum becomes sharper due to the 
weakening of quasiparticle scattering near the normal- 
metal/superconductor micro-constriction and the barrier 
strength also decreases. Further pressure of the tip on 
the sample surface may simply flatten the point (vali- 
dating more real point contact channels) over the same 
minimal thickness of a tenacious barrier layer. In this 
case, if no any artificial effect shows up, there should 
be no obvious change on the measured spectrum in a 
range of junction resistance (i.e., the normalized spectra 
should be identical) until the junction is damaged even- 
tually. All measurements are carefully checked to be in 
such a regime. We have also found that the determined 
gap value is independent on the barrier strength and the 



measured locations (taking into account the slight varia- 
tion of local T c ), indicating that we stood a good chance 
to detect the bulk properties of the samples. 

The temperature dependence of superconducting gap 
are presented in Fig. [3]j) in a reduced scale. The uni- 
versal A(T) relation was found for all studied doping 
levels in good agreement with the BCS theory. The 
derived A(0) and T c for various locations on different 
samples are summarized in Fig. [4] For comparison, we 
also replot the data in previous reports (grey symbols) 
[S d) EH EE E3 for optimally doped samples and an 
overdoped one. The coupling ratio is almost a constant 
with a value around 1.7. If the non-monotonic d-wave 
model is used, this ratio becomes about 1.6 (as exem- 
plified by the open circles). Both cases belong to the 
weak coupling regime. As an example to demonstrate 
the validity of the methodology for determining the su- 
perconducting gap, we also presented in Fig. 2] the data 
from Nb-tip/Au-foil point contacts using the same ex- 
perimental apparatus. The determined coupling ratio is 
consistent with the previous reports [29[ and BCS theory 
with a high precision. 



IV. SUMMARY 

In summary, by investigating the point contact spectra 
of the electron-doped cuprates, we show that the temper- 
ature dependence of the superconducting gap follows the 
BCS prediction very well in a wide doping regime with a 
universal weak coupling ratio of A(0)/fcsT c = 1.72±0.15. 
Therefore, the electron-doped cuprates are weak cou- 
pling BCS superconductors although the non-monotonic 
d-wave pairing symmetry may be favorable. 

Note added: During the preparation of this 
manuscript, we became aware that a recent report in- 
dicates the weak coupling BCS dirty superconductivity 
in an electron doped cuprate based on the measurements 
of SIS' junctions[34j. We also became aware of the recent 
STM paper (to be published) reporting larger gaps than 
that in this work. 
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